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Abstract 


Recent developments from a continuing effort 
to provide an equivalent plate representation for 
aircraft structural analysis are described. 
Previous work provided an equivalent plate 
analysis formulation that is capable of modeling 
aircraft wing structures with a general planform 
such as cranked wing boxes. However, the 
modeling is restricted to representing wing boxes 
having symmetric cross sections. Further 
developments, which are described in this paper, 
allow modeling of wing cross sections having 
asymmetries that can arise from airfoil camber or 
from thicknesses being different in the upper and 
lower cover skins. An implementation of thermal 
loadings, which are described as temperature 
distributions over the planform of the cover 
skins, has been included. Spring supports have 
been added to provide for a more general set of 
boundary conditions. Numerical results are 
presented to assess the effect of wing camber on 
the static and dynamic response of an example 
wing structure under pressure and thermal 
loading. These results are compared with results 
'rom a finite element analysis program to 
indicate how well a cambered wing box can be 
represented with an equivalent plate formulation. 


Nomenclature 


a,b,c,e,f ,g 

A i’Y C k 


h 

K 


ij 


P ,P ,P 
u V w 



t 

T 

U,v,w 


x.y 


V*i 


-Planform dimensions (see Fig. 3) 
-Coefficients of polynomial 

displacement functions for U, V, 
and W components, respectively 
-Wing box depth 
-Stiffness submatrix for i,j 

displacement components 
-Load vectors corresponding to U,V,W 

deformations 

-Lamina stiffness matrix for k th 
layer 

-Components of lamina compliance 

matrix (see eqn.(6)) 

-Thickness of cover skin layer 
-Temperature 

-Displacement functions in the x,y,z 
directions, respectively 
-Global streamwise and spanwise 
coordinates, respectively 
-Polynomials in x and y for defining 

displacement functions 
-Location of mid-camber surface 
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-Thermal coefficients of expansion 

in x-y coordinates 
-Thermal coefficients of expansion 

along orthotropic material axes 
-Total strains in x-y coordinates 

-Stress induced strains along 

orthotropic material axes 
-Temperature induced strains along 

orthotropic material axes 
-Local nondimensional streamwise and 
spanwise coordinates, respectively 
-Function used in surface fitting 
procedure (see eqn.(15)) 

-Potential energy 
-Stresses in x-y coordinates 

-Orientation of orthotropic material 
axes (see Fig. 1b) 

-Parameter to define magnitude of 
wing camber 


Introduction 


Simplified beam or plate models of aircraft 
wing structures are often used for analysis 

during early preliminary design. 1 For example, 
an equivalent plate model of the wing structure 
is used in the TSO (Aeroelastic Tailoring and 

2 3 

Structural Optimization) computer program. 

This program has had widespread use for 

4 

aeroelastic tailoring of composite wings. 

However, the structural analysis formulation used 
in TSO is limited to trapezoidal planforms. 

Recently, a new equivalent plate analysis 
formulation has been developed with the 
capability to model aircraft composite wing 
structures with general planform geometry such as 

cranked wing boxes. ^ This formulation contains 
only bending terms in the expression of the 
equivalent plate energy. Therefore, the modeling 
is restricted to representing wing boxes having 
cross sections that are symmetric about a 
midplane. Also, all applied loads, such as 
aerodynamic pressures , must act normal to this 
midplane. 

The present paper describes a further 
generalization of the equivalent plate 
representation. This extension allows modeling 
of unsymmetric wing cross sections which can 
arise from airfoil camber or from having 
different thicknesses in the upper and lower 
cover skins. This capability is included in the 
formulation by adding stretching terms to the 
deformation of the reference surface which does 
not have to be located at the midplane of the 
wing. The addition of stretching deformations in 
the reference surface provides fully-coupled 
bending-stretching behavior of a composite plate 
to be represented. 
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Inplane loads can also be applied in this 
generalized formulation. An implementation of 
thermal loadings, which are described as 
temperature distributions over the planform of 
the cover skins, has been Included. Spring 
supports have been added to provide for a more 
general set of boundary conditions. 

The paper will include a description of the 
analytical formulation along with an example to 
demonstrate the additional features. Results 
will be presented to assess the effect of wing 
camber on the static and dynamic response of an 
example wing structure under pressure and thermal 
loading. These results will be compared with 
corresponding results from a finite element 
analysis program to indicate how well a cambered 
wing box can be represented with a flat plate 
formulation. 

It is assumed that the reader of this paper 
is familiar with reference 5 which contains a 
description of the basic method. The material 
contained in the present paper will focus on 
providing a description of extensions to the 
previous work and will repeat only a minimal 
amount of information contained in reference 5. 


Analytical Modeling 


The wing box structure is represented as an 
equivalent plate in this formulation. Planform 
geometry of this equivalent plate is defined by 
multiple trapezoidal segments as Illustrated by 
the two-segment box in Figure la. The global 
cartesian coordinate system has the x-axis in the 
streamwise direction and the y-axis in the 
spanwlse direction. A separate local coordinate 
system is associated with each segment. These 
local coordinates are nondimens lonallzed such 
that £ refers to a fraction of the local chord 
and n refers to a fraction of the span for a 
given segment as indicated in Figure 1b. The 
subscripts on the £ and n coordinates, shown in 
Figure 1 to refer to a particular segment, are 
omitted in the remainder of this paper since the 
development of the analysis method is described 
for a typical segment. The orthotropic material 
axes, denoted 1 and 2, are also shown in Figure 
1b for a typical layer in the cover skin. 

The cross-sectional view of a typical 
segment in Figure 2 Illustrates the analytical 
modeling of the wing box structure. The location 
of the mid-camber surface of the wing Is defined 
as the distance, zc, from a user specified 
reference plane. This distance varies over the 
planform of each segment and is expressed as a 
polynomial in the global coordinates x and y. 


zc(x,y) 


W z 10 x+z 20 x 


m n 

.♦z x y 
mn 


( 1 ) 


The coefficients z are constants which are 
mn 

defined by the analyst f or each segment. 
Similarly, the depth of the structural box, which 
also varies over the planform, is defined by the 
analyst again in polynomial form. 


with the thickness of each layer being defined 
independently also in polynomial form 

V X>y) * t 00 H 10 Xtt 20 x2+t 0l jrt -V" yn (3) 

The properties of the layers can be defined to 
represent wing skins which are stiffened panels 
or composite laminates. Orientation of the 
stiffness properties and corresponding thickness, 
are specified for each layer and the orientations 
and thicknesses can be different in different 
planform segments. In the present 
implementation, the thickness of a layer in the 
lower skin is specified as a factor times the 
thickness of the corresponding layer in the upper 
skin. This simplification is implemented to 
reduce the number of variables which must be 
considered in a design application of this 
procedure. The exponents of the polynomial 
terms, m and n, are not necessarily the same for 
equations (1), (2), and (3). These exponents are 
specified by the analyst and values of the 4th 
power have been found sufficient for representing 
a fairly general class of wing box geometry. 

New options have been added for modeling the 
volume between the covers. This volume can be 
(a) considered empty to represent a dry wing, (b) 
used to represent fuel by assigning a fuel 
density, or (c) specified as solid structure by 
giving material stiffness properties. Option 
(c), is provided to analyze wind tunnel models 
for conditions where aeroelastic effects may be 
important. Although not shown in Figure 2, rib 
and spar caps can be represented as axial members 
and concentrated masses can be defined for 
vibration analysis as described in reference 5. 

Thermal loading is included for static 
analysis along with concentrated forces and 
distributed pressure loads. Thermal loads are 
specified as temperature distributions over the 
planform of the cover skins. The temperature is 
assumed constant through the thicknesses of all 
layers making up the skin. In the present 
procedure, material properties in the skins are 
not considered to be a function of temperature. 


Energy Expression for Plates 


An expression f or the strain energy of a 
plate segment is needed for the Ritz solution 
procedure that is used. This expression is 
developed in terms of bending and stretching of 
the reference plane shown in Figure 2. The 
Klrchhoff assumption is made that lines normal to 
the reference plane remain straight and normal 
under deformations of the equivalent plate. This 
condition is analogous to the assumption that 
cross sections of a beam remain plane and normal 
to the centroidal axis in beam theory. The 
deformations in the equivalent plate are given as 


U - U n - z W, 

0 x 


V n - z W, 

0 *y 


(4) 


hU.y) - \ o ^ lo x^ 0 x 2 *h 0I y*...*h m *V 
The cover skins consist of orthotropic layers 


( 2 ) 


W - W 

where U Q and V Q refer to stretching of the 
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Analyslg Procedure 


reference surface in the x and y directions at 
z*0, and W is a deformation normal to the surface 
caused by bending. Corresponding strains in the 
x and y directions are given by 

e - +0 - z W, 

x Crx xx 


+ Vy 


Vy + Vx 


2 z W, 


In order to include thermal loading, the 
total strain, e, is separated into the sum of 
stress-induced components, e 1 , and temperature- 
induced components, e" as in reference 6. These 
strains can be expressed for each layer in the 
cover skins in terms of orthotropic material 
properties along the axes 1 and 2 as 


S 1 2 S 22 


where the first and second set of terms In the 
summation are e f and e" respectively. The 
stresses and stress-induced strains in each 
lamina are transformed to the global (x,y) 
coordinate system and used in the expression for 
strain energy as 


n - 1/2 / _ (o e ’ 

vol x x 


o e ' ♦ o e 1 ) dV 
Y Y xy xy 


This energy expression can be written in terms of 
total strains as the 3um of two parts n - ♦ 


n, - 1/2 ; vol ( £ } T C Q k ] { e ! dV 


-1/2 / yol { a } T [ Q k ] { e ) dV 


T 

where { e } » { e e e } are total strains, 

x y xy 

T 

{ a } - { a a a } are ”ef fective" 

x y xy 

thermal coefficients as discussed in reference 7, 

k 

and [ Q ] is a lamina stiffness matrix. 

Equation (8a) gives the strain energy of the 
plate and equation (8b) Is the potential energy 
of the thermal loading. When the strain energy 
of the plate from equation (8a) is expanded in 
terms of the reference surface displacements in 
equation (5), a total of 29 terms result as given 
in reference 8. Note that for bending alone 
there are only 6 terms in the energy expression 
as given in reference 5. Therefore, the addition 
of the stretching deformations to handle cambered 
wing boxes requires considerable additional 
computation. 


The Ritz method is used to obtain an 
approximately stationary solution to the 
variational condition on the energy of the 
structure and applied loads. In this application 
of the Ritz approach, each component of the 
deformation of the reference surface is assumed 
to be the sum of contributions from sets of 
specified displacement functions 

U - f A 1 X(x) i Y(y) 1 

V - jl B j X(x) J Y(y) j (9) 

w - $ C k X(x) k Y(y) k 

As in reference 5, these displacement functions 
are specified as products of terms from a power 
series in the x-direction with terms from a power 
series in the y-direction. Different sets of 
functions can be used for the U, V, and W 
components of the deformation. 

For static analysis, the Ritz procedure as 
described In reference 5 produces a system of 
simultaneous equations which can be solved for 
the unknown coefficients in equation (9) to 
minimize the total energy expression. These 
equations can be written in partitioned matrix 
form as 



These equations represent the fully-coupled 
bending-stretching behavior of a composite plate. 

The evaluation of the terms in the 
submatrices of the stiffness matrix of equation 
(10) requires the integration of the strain 
energy in equation (8a) over the volume of 
structural material. This integration involves 
functions of the camber shape, zc, the wing 
depth, h, and the thickness of each cover skin 
layer, t k , along with the corresponding lamina 

k 

stiffness matrix, Q . To simplify calculations 
and make the resulting stiffness matrix a linear 
function of layer thickness, all layers in a skin 
are assumed to be located at the same distance 
from the reference surface. During integration 
their thicknesses are simply summed to yield the 
appropriate stiffness values. This 
simplification is invalid for plates with a small 
distance between the covers but results in only a 
small error for typical wing structures where the 
depth of the wing is large compared to the 
thickness of the covers. Exact integral 
expressions are formed and evaluated at limits on 
the upper and lower covers for the material 
between the covers. The expressions describing 
the model shown in equations (1)-(3) must be 
combined with quadratic functions of the strains, 
hence displacement functions, as given in 
equation (5) to complete the expression for 
strain energy of the plate as shown in equation 
(8a). The formation and evaluation of the terms 
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In these integral expressions is algebraically 
cumbersome and involves tedious manipulations. 
However, the mathematical operations involved 
yield polynomials in x and y and the same library 
of subroutines which were developed in reference 
5 are used to perform all the mathematical 
operations on these polynomials in an efficient 
manner. 


Thermal Loads 

Thermal loading is included for static 
analysis along with concentrated forces and 
distributed pressure loads. Thermal loads are 
specified as temperature distributions over the 
planform of the cover skins. The temperature is 
assumed constant through the thicknesses of all 
layers in the skin. The thermal loads are 
calculated using equation (8b). The thermal 
coefficients which are needed for each layer are 
transformed to the global x-y coordinate system 
by 

2 2 

a x * cos e a 1 + sin 0 a 2 

2 2 

dy - sin e a j + cos 0 (11) 

a xy - 2 sine cose (o^-dg) 

where 0 is the angle from the x-axls to the 1- 
axis of the orthotropic material. At present, 
the stiffness properties of the layers, Q, are 
not considered to be a function of temperature in 
the analysis procedure. 

If the temperature distribution, AT(x,y), is 
input as a polynomial, the necessary integrations 
are performed using closed form expressions. 
Often, the temperatures have been calculated in a 
separate program on a grid over the wing 
planform. In this case, numerical Integration 
can be performed with the appropriate values of 
temperature, thermal and stiffness coefficients, 
strain, and planform area associated with each 
point in the grid being used. This integration 
must be performed for each of the displacement 
functions which are used in the analysis as 
described in reference 5. This procedure is 
analogous to the procedure used for application 
of aerodynamic pressure loads in the equivalent 
9 

plate formulation. 

The capability to consider both aerodynamic 
and thermal loading Indicates a potential use of 
this procedure for static aerothermoelastic 
calculations for advanced vehicles where 
aerodynamic heating is an Important design 
consideration. Such calculations would involve 
an iteration on a combination of thermal loads, 
structural deformations, and aerodynamic loads to 
arrive at a set of loads which are consistent 
with the deflected shape of the loaded wing at 
some specified flight condition. 


Definition of Constraints 

For static analysis, rigid body motion of 
the plate must be constrained. These constraints 
are often referred to as boundary conditions. In 


reference 5, the example equivalent plate wing 
model was clamped along the x-axis which 
represented an aircraft centerline. This is 
accomplished by excluding selected terms from the 
set of displacement functions as follows: 

(a) the condition that the displacement, W is 
zero at y-0 can be specified by excluding 

all y° terms from the set of displacement 
functions 

(b) the condition that the slope, W, y is zero 

at y-0 can be specified by excluding all 

y terms from the set of displacement 
functions 

(c) clamped conditions at y-0 is specified by 
imposing both (a) and (b). 

This method of specifying constraints is limited 
to the x and y axes. Other methods of applying 
constraints which are available in the 

literature were evaluated as candidates to 
provide for specification of more general types 
of boundary conditions. The Lagrange multiplier 
method allows exact specification of constraints 
but at the expense of introducing additional 
equations with zeros on the diagonal. The method 
selected for use is an adaptation of the penalty 
function method discussed in reference 10. This 
method will satisfy the constraints only 
approximately. A penalty number is used and the 
larger the penalty number the better the 
constraints will be achieved. For the purposes 
of constraining the equivalent plate, this 
penalty number is taken to be the stiffness of a 
spring at the location of the desired constraint. 
Such springs can resist translation or rotation 
and are defined at a point or distributed along a 
specified line. Although displacements cannot be 
specified to be exactly zero at a selected 
location, use of sufficiently stiff springs will 
provide a good approximation to the desired 
condition. 

A combination of excluding terms from the 
displacement functions and application of stiff 
springs is used to define the constraints on the 
example wing which will be discussed in a later 
section of this paper. The input data used to 
define the displacement functions for the 
analysis are 


DIR 

NX 

X powers 


NY 


Y 

powers 



1 

5 

0 1 2 

3 

T 

”T> 

0 

1 

2 

3 

~r 

5 

2 

5 

0 1 2 
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5 

1 

2 

3 

k 

5 


3 

5 

0 1 2 

3 

k 

6 

0 

2 

3 

4 

5 
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The directions 1,2,3 refer to the U,V,W 
deformations and NX, NY indicate the number of 
terms to be used in the x and y directions 
respectively. The displacement function in each 
direction is composed of a polynomial containing 
the sum of terms from all combinations of the x- 
powers and y-powers which are given. Constraints 

for the example are imposed by omitting y° for 
the U-deformation to provide a deflection 

constraint at y-0 and by omitting y 1 for the W~ 
deformation to constrain the slope at y-0. 
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In addition to the more general capability 
for specifying constraints, an elgenproblem shift 

parameter 11 has been Included for use In 
vibration analysis. This shift parameter allows 
a vibration' analysis to be performed on a model 
with unconstrained (rigid body) motions. 


Implementation of Method 

This equivalent plate analysis procedure is 
being developed for use with mathematical 
optimization procedures for application in early 
preliminary design. Therefore, an important 
facet of the development is the implementation 
into a computationally efficient computer 
program. Although a detailed description of the 
computer program is beyond the scope of this 
paper, some of the significant approaches used to 
achieve the desired efficiency are outlined in 
this section. 

Generation of Stiffness Matrix 

The terms associated with calculating 
coefficients of a stiffness matrix for an 
anisotropic plate segment are algebraically 
cumbersome and tedious to manipulate. The 
general procedure for calculating these terms for 
only the bending deformation, W, is described in 
reference 5. To include the stretching 
deformations U and V, the size of the stiffness 
matrix is increased as indicated in equation 
(10). The procedure involves evaluation of the 
integral shown in equation (8a). The strains 
needed for this evaluation can be written as the 
sum of components from the U, V, and W 
deformations as shown in equation (4). Each 
submatrix in equation (10) is formed by using the 
appropriate combination of the components in the 
various directions. 

The procedure is Implemented by generalizing 
the corresponding procedures used for in 

reference 5. These calculations can then be 
performed by looping through this common set of 
generalized procedures with the origin of each 
submatrix and the definition of the strain 
components being changed at the beginning of each 
loop to correspond to the desired combination of 
the U, V, and W directions. The strain 
components are defined in a table containing the 
direction, the coefficients, and an indicator to 
specify the proper derivative operation. 

Integral Tables 

The evaluations outlined above are performed 
for all combinations of terms which are specified 
in the displacement functions for each direction 
in order to complete a submatrix. All these 
evaluations involve generation and integration of 
lengthy polynomial expressions. Since these 
expressions are all sums of terms containing a 

coefficient multiplied by x m y n , tables are 
generated with each entry containing the Integral 
of such terms over the planform of the segment 
being evaluated for all combinations of m and n 
that are required. Entries in the table are 
evaluated by integrating over the unit square for 
each plate segment after expressing the terms in 


the local coordinate system. The coordinate 
transformations are given in terms of the 
planform variables shown in Figure 3 as 

x - e + a£ + (f-e ) rj + (c-aKn (12) 

y - g + bn (13) 

The coordinates of the differential area are 
transformed by the determinant of the Jacobian as 

dxdy - [ ab + (c-a)bn ] d£dn (14) 

The distance from the reference surface to 
either the upper or lower cover skin is indicated 
by the symbol z in the expression for the W 
component of strain in equation (5). This 
distance can be expressed in polynomial form in 
terms of the mid-camber and wing depth 
definitions given in equations (1) and (2). The 
maximum power of z is 2 for the pure bending, K ww 

submatrix for cover skins and is 3 after 
integrating through the wing depth for solid 
plates. In the implementation used herein, eight 
integral tables are generated corresponding to 
the 0, 1, 2, and 3 powers of z for the upper and 
lower surfaces. The tables contain evaluation of 

terms corresponding to z^ times x m y n which are 
used to generate the stiffness submatrices for 
appropriate combinations of strain components. 

Input of Analytical Model 

The definitions of the mid-camber surface, 
wing depth, and thicknesses of layers in the 
cover skins are given in terms of the global 
coordinates in equations (1)-(3). Such a 
definition results in all subsequent polynomials 

that are generated being the sum of x m y n terms. 
Retaining this form allows the efficient handling 
of the integral operations as described above. 

An input option is provided for defining the 
model directly in the x-y coordinate system. 
However, it is often more convenient for the user 
to input these quantities in the local £-n 
coordinate system of each segment as was done in 
reference 5. Thi3 input option has been included 
in the present implementation but a 
transformation from the £-n system to the global 
x-y system is necessary. A direct transformat ton 
using the inverse of equations (12) - (14) 
results in a complicated form rather than the 
desired simple power series polynomials. 
Therefore, an approximate transformation to the 
power series polynomials is used. 

This approximate method minimizes the square 
of the difference between the input polynomial in 
the system and a power series polynomial in 
the x-y system integrated over the planform of 
the segment 1 . The function to be minimized can be 
written as 


■ 2 •'area* 1/2 ° 2 * FG + 1/2 p2 > dA (15) 

where F - f(£»n) which is input 

and G - g(C mn x m y n ) is the desired polynomial. 
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Differentiating equation (15) with respect to the 

unknown coefficients, C , yields a set of 
mn 

simultaneous equations which can be generated and 
solved using the same procedures at that for the 
potential energy given in equation (8). This 
solution is the desired surface fit and the 
procedure can be used to obtain the coefficients 
for equations (1)-(3) to define an analytical 
model . 

Analysis Time 

The time required to perform an analysis 
includes the time required for input preparation 
to define the analytical model, the execution 
time of the program, and time required for the 
analyst to assimilate the output data. The 
specification of model characteristics as 
continuous distributions in polynomial form 
requires only a small fraction of the volume of 
input data for a corresponding finite element 
structural model where geometry and stiffness 
properties are specified at discrete locations. 
Equivalent plate models are applicable for cases 
where the stiffness can be represented as a 
continuous distribution within segments, 
therefore limiting the complexity of modeling. 

The execution time Is a function of 
complexity and refinement of the analytical 
model. The number of degrees of freedom used in 
an analysis to give adequate accuracy is problem 
dependent. Results are presented in this paper 
for two different sets of displacement functions 
for the equivalent plate model and for two 
different element mesh refinements for the finite 
element model. Computational times and selected 
results are compared to give an indication of the 
accuracy and computational efficiency of these 
procedures for different levels of modeling. 

The number of degrees of freedom in an 
equivalent plate model is readily changed by 
simply changing the specification of terms in the 
displacement functions rather than refining the 
number of joints and elements in a finite element 
model. The output quantities such as 
displacements, strains, stresses, and vibration 
mode shapes are calculated on a user specified 
grid over each plate segment. Since the 
definitions of these quantities are given as 
continuous polynomials, the grid can be as fine 
or as coarse as desired. Presently, no graphical 
presentation of output data has been Implemented. 


Application and Results 
Model Description 

The planform of an example wing box which is 
being used to evaluate the effects of including 
camber in the formulation is shown in Figure 4 . 

A two segment planform, representative of a 
fighter aircraft wing, is used for testing 
purposes. The planform is composed of a clipped 
delta outer segment with a 45 degree leading edge 
sweep and an inner segment to represent a carry 
through structure. A delta planform is selected 
for evaluation because cambering of the long 
chord length in the inboard portion of the wing 
leads to significant changes in the z coordinate 
of the wing surfaces. Two loading conditions 


were used for testing. The first is an 
externally applied loading of a uniform 1 psi 
pressure acting over the planform of the outer 
segment. The second is a thermal loading which 
is specified as the temperature of the lower 
surface of the wing segment being 100 degrees F 
greater than that of the upper surface. Boundary 
conditions are Imposed by: (a) by excluding 
selected terms from the displacement functions as 
described in an earlier section to constrain the^ 
slope of the bending deflection, W, and the 
spanwise stretching deflection to zero at y-0 qnd.'* 
(b) using stiff springs to constrain the U and W 
deflections at support point A and the W 
deflection at. support point B as shown in Figure 
4. 

A simple expression for the shape of the 
mid-camber surface is used to evaluate the 
effects of camber. This shape is specified as a 
quadratic function of the chord length (parabolic 

arc) with the maximum dimension, zc , located 

max 

at 50% chord location. The amount of camber will 

be defined as $; where $ - zc /local chord. 

max 

Results are presented for no camber (symmetric 
wing) and a camber of $-.03. The equivalent 
plate analysis for the cambered wing is performed 
using two different sets of displacement 
functions for the bending deflection, W. The 
first set, referred to as the 4x6 analysis, has 
the exponents shown earlier in the discussion of 
constraints with 4 being the largest exponent in 
the x-direction and 6 being the largest exponent 
in the y-direction. In the second set which is 
referred to as the 5x5 analysis, 5 is the largest 
exponent in both the x- and y-directions . 

Results from the equivalent plate analysis 
for the cambered wing are compared with 
corresponding results from the EAL finite element 

analysis program. The EAL model is built up of 
membrane rib, spar, and cover elements with the 
grid of cover elements shown in Figure 4. 

Analyses are performed for two levels of 
modeling. One model has a grid of cover elements 
with twice the refinement as the other model in 
both the x and y directions and has additional 
rib and spar elements corresponding to the edges 
of each cover element. The models will be 
referred to as the 435 degree of freedom (d.o.f.) 
model and the 1565 degree of freedom model as 
indicated in Figure 4. 

Numerical Results 

An Indication of the relative level of 
accuracy of results and corresponding 
computational time for the different models used 
to analyze the cambered wing box structure are 
presented in Table 1 . A measure of the overall 
accuracy is given by the total strain energy of 
the deformed structure for the pressure and 
thermal load cases. Increased accuracy is 
indicated by an increasing value of strain 
energy, so that the 1565 degree of freedom finite 
element model is shown to provide the best 
overall results in this study. The corresponding 
increase in computational time required to 
provide the improved accuracy is also shown in 
Table 1 . The equivalent plate analyses required 
less computational time than the finite element 
analyses. However, the magnitude of the 
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difference is very dependent on the level of 
modeling used. Values of deflections and 
stresses at selected locations on the wing along 
with vibration frequencies are presented in the 
remainder of this section to provide a detailed 
comparison of results. 

Static displacements for the uniform 
pressure load are presented along the trailing 
edge of the wing in Figure 5. These results 
indicate that camber reduces the deflection by 
approximately 5 % at the wing tip. This 
relatively small difference suggests that a flat 
representation of a cambered wing may be adequate 
for aeroelastic calculations. The difference 
between the deflections from the 4x6 and the 5x5 
equivalent plate models are negligible so that 
only results from the 4x6 model are shown in 
Figure 5. The deflection obtained from the 
finite element model of the cambered wing with 
435 degrees of freedom is in close agreement with 
that obtained from the equivalent plate method. 
However, these deflections are lower than those 
given by the 1565 degree of freedom finite 
element model. 

Distributions of the stress in the spanwlse 
direction for the upper surface and lower surface 
of the wing are shown in Figures 6-8. These 
distributions are given at the y * 20, 80, and 
120 locations illustrated on Figure 4. In some 
regions there are significant differences between 
the stresses for a flat wing and the stresses for 
a cambered wing. For example, this difference is 
as large as a factor of 2 at some points of the 
y-120 location. These results Indicate the 
importance of considering camber if the analysis 
is used in a procedure to size cover skins for 
strength requirements since the resized 
thicknesses are proportional to the calculated 
stresses. 

In general, the agreement in stresses for 
the equivalent plate procedure and the finite 
element analysis Is acceptable for use during 
preliminary design. The largest difference 
occurs in the bottom surface at the y»80 
location. This is the region of transition from 
•the carry through structure to the outer delta 
portion and It appears that the degree of the 
polynomials used for the displacement functions 
is not large enough to represent the gradients 
occuring in this region. The 5x5 equivalent 
plate model gives better agreement than the 4x6 
model in this region. Attempts to use larger 
exponents for this particular problem gave ill- 
conditioned equations which caused the solution 
subroutine to abort with an error message. 

The deflections along the leading and 
trailing edges of the wing box produced by a 
temperature differential of 100 degrees F on the 
lower surface of the wing are shown in Figure 9. 

A negative 0.6 degree angle of attack is produced 
at the wing tip. Such twisting of the wing that 
occurs from this thermal loading could affect 
aerodynamic loads and indicates a potential use 
of this procedure for aerothermoelastic 
calculations. 

Vibration frequencies calculated for the 
cambered and uncambered wings are shown in Table 
2. Camber was found to have only a small effect 
on the values of the first ten frequencies. The 
values of frequencies from the equivalent plate 
analyses are higher than from the finite element 
analyses, indicating that the displacement 


functions which were used did not represent the 
vibrational behavior of the wing as well as the 
finite element model. However, the level of 
accuracy of the frequencies indicate that the 
equivalent plate method provides an acceptable 
representation of the vibration characteristics 
of the wing for use during preliminary design. 


Concluding Remarks 

A description is given of a further 
generalization of the equivalent plate 
formulation to provide capability to model 
aircraft wing structures with unsymmetric cross 
sections. The analytical procedures used to 
provide this capability are given along with some 
of the methods used for implementing these 
procedures into a computer program. 

Some typical numerical results are presented 
to assess how well a cambered wing box can be 
represented by a plate formulation. In general, 
the difference in results between static 
displacements and vibration frequencies for the 
cambered wing example and a flat representation 
is less than 1 0% • However, the difference in 
stresses is significant; a factor of 2 difference 
in some regions. This large difference in 
stresses indicates the importance of including 
the effects of camber in a general structural 
resizing procedure. 

The degree of the polynomials used for the 
displacement functions was limited by ill- 
conditioning of the resulting equations and the 
same level of accuracy of the finite element 
results could not be achieved with the equivalent 
plate procedure. However, the results from the 
equivalent plate analyses exhibited the same 
trends and approached the accuracy of the finite 
element analysis results. Considerably less 
total analysis time was required to generate the 
equivalent plate results. 

In summary, the application o'* the 
generalized equivalent plate analysis procedure 
to cambered wing structures is shown to produce 
results with adequate levels of accuracy In a 
shortened analysis time. Hence, the procedure 
provides a useful capability r or the analysis of 
aircraft structures during early preliminary 
design. 
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Table 1 . Comparison of strain energy and computational time 
for different levels of structural modeling. 


Model 

Strain Energy, 
Pressure 
Load 

in. -lb. 
Thermal 
Load 

Computation 
Time, sec. 

Eq. plate 
4x6 

19660.7 

7887*1.0 

50.9 

Eq. plate 
5x5 

19700.1 

78917.6 

5*1.7 

FEM 

435 dof 

19919.9 

8121 1 .5 

122.0 

FEM 

1565 dof 

21216.6 

81897.9 

588.0 


Table 2. Comparison of natural frequencies 
from vibration analysis. 


Number 

Eq. plate, 4>«Q 
4x6 model 

Eq. plate 

. *-.03 

FEM. 

4>-.03 

4x6 model 

5x5 model 

435 do' - 

1565 dof 

1 

1 6.2*4 

16.52 

16.52 

16.38 

15.90 

2 

53. *16 

54.54 

54.41 

51 .89 

52.58 

3 

80.63 

77.23 

77.17 

72.27 

72.56 

4 

83.01 

82.01 

82.06 

79.63 

81.08 

5 

88.65 

93.72 

93.12 

80.58 

82.20 

6 

127.54 

128.23 

128.83 

119.43 

124.28 

7 

184.06 

185.06 

180.68 

173.09 

173-30 

8 

192.34 

193.81 

191.36 

178.77 

185.71 

9 

212.40 

213.60 

207.06 

199.33 

199.57 

10 

227.48 

231.32 

233.87 

218.25 

214.80 


8 





Eq. plate 4 x 6, <t> = 0 
Eq. plate 4x6, 

Eq. plate 5x5 . 

FEM, 435 d.o.f. 9 ' 

FEM, 1565 d.o.f. 


0= .03 


Eq. plate 4 x 6, <1> = 0 
Eq. plate 4 x 6, \ 

Eq. plate 5 x 5 ( 
FEM, 435 d.o.f. 9 
FEM, 1565 d.o.f.) 


Figure 6a. Stress distribution at y-20; 
top cover. 


Figure 6b. Stress distribution at y-20; 
bottom cover. 


Eq. plate 4x6, 4>=0 

Eq. plate 4x6 

Eq. plate 5x5 m 

FEM. 435 d.o.f. 

FEM, 1565 d.o.f. 


A 

<D= .03 O 

J 


Eq. plate 4 x 6, d) = 0 

Eq. plate 4x6 

Eq. plate 5x5 

O FEM, 435 d.o.f. 

A FEM, 1565 d.o.f. 

A ^-4-. , 


4)= .03 


Figure 7a. Stress distribution at y-80; 
top cover. 


Figure 7b. Stress distribution at y-80; 
bottom cover. 





Eq. plate 4 x 6, <t> = 0 



$ 
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Figure 8a. Stress distribution at y-120; 
top cover. 


Figure 8b. Stress distribution at y-120; 
bottom cover. 



Figure 9. Thermal deflection of wing box. 
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